Bounds on neutron-star moments of inertia and the evidence for general relativistic 

frame dragging 
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Recent X-ray variability observations of accreting neutron stars may provide the first evidence for 
frame dragging effects around spinning relativistic objects. Motivated by this possibility and its 
implications for neutron-star structural properties, we calculate new optimal bounds on the masses, 
radii, and moments of inertia of slowly rotating neutron stars that show kilohertz quasi-periodic 
oscillations (QPOs). These bounds are derived under minimal assumptions about the properties of 
matter at high densities and therefore are largely independent of the unknown equation of state. 
We further derive a semi-analytical upper bound on the neutron-star moment of inertia without 
making any assumptions about the equation of state of matter at any density. We use this upper 
bound to show that the maximum possible nodal precession frequency of an inclined circular orbit 
around a slowly spinning neutron star is i/np — 45.2(^/300 Hz) Hz, where u B is the spin frequency 
of the neutron star. We conclude that the nodal-precession interpretation of low-frequency QPOs 
in accreting neutron-stars is inconsistent with their inferred spin frequencies or the identification of 
the highest-frequency QPO with that of a circular Keplerian orbit in the accretion disk. 
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I. INTRODUCTION 

Astrophysical observations of general relativistic phe- 
nomena are rare and often confined to effects observable 
within our solar system Q. Among galactic sources, rela- 
tivistic compact objects, such as neutron stars and black 
holes, are the prime candidates for the detection of such 
phenomena. In these sources, however, the effects of gen- 
eral relativity are often convolved in complex ways with 
other astrophysical phenomena making their identifica- 
tion often impossible. The observation, via precise pulsar 
timing, of orbital shrinkage in double neutron-star sys- 
tems represents the most successful to date identification 
of pure general relativistic effects in galactic compact ob- 
jects §. 

The recent discovery with the Rossi X-ray Timing Ex- 
plorer of quasi-periodic oscillations with variable kilo- 
hertz frequencies (hereafter kHz QPOs) from many ac- 
creting neutron stars has provided us with a new clock 
that measures accurately their variability The cen- 
troid frequencies of the highest-frequency kHz QPOs in 
each source have been identified with Keplerian frequen- 
cies of stable orbits in the accretion disks very close to 
the stars These QPOs can therefore be used as 

probes of the physical conditions in regions where gen- 
eral relativistic effects in the strong-field regime are non- 
negligible l|. This identification has lead to new upper 
bounds on the masses (typically < 2.2 M Q ) and radii 
(typically < 15 km) of the accreting neutron stars p,||. 
In one specific X-ray source (4U 1820—30), the associa- 
tion of a kHz QPO with the Keplerian orbital frequency 
at the radius of the innermost stable circular orbit around 
the neutron star may provide us with the first evidence 
of this prediction of general relativity as well as with the 



identification of the first relatively massive (~ 2.2 Mq) 
neutron star 0. 

In many accreting neutron stars, a third low-frequency 
(~ 10 — 70 Hz) QPO, the so-called horizontal-branch os- 
cillation (HBO), is often observed simultaneously with 
the kHz QPOs [||||. Soon after its discovery the ori- 
gin of this QPO was linked to the interaction between 
the accretion disk and the magnetosphere of the neu- 
tron star Recently, however, this interpretation has 
been challenged |l(| by the understanding that the accre- 
tion disks appear to penetrate closer to the stellar surface 
than what is required by the model H^JnJ. As an al- 
ternative model, Stella & Vietri |l(J suggested that the 
frequency of the HBO is the general relativistic Lense- 
Thirring precession frequency, caused by frame dragging, 
of an inclined circular orbit that has a Keplerian orbital 
frequency equal to the frequency of the highest-frequency 
kHz QPO. The physical mechanism responsible for pro- 
ducing brightness oscillations at the precession frequency 
of a particular orbit in the accretion disk is still a matter 
of active research |l^]. However, when the effects on the 
precession frequency of the quadrupole moment of the 
stellar gravitation field can be neglected, the theoreti- 
cally predicted correlation between the precession and 
orbital frequencies is consistent with observations [pTo JTT | 
(see Ref. |llj for a discussion of the discrepancy between 
the predicted and observed trend at high frequencies). 
The normalization of the observed correlation, though, is 
relatively high. In other words, the magnitude of the ob- 
served HBO frequencies can be accounted by the model, 
only if the moments of inertia of the neutron stars are 
~ 4 — 5 times larger than predicted by any realistic equa- 
tion of state [jn|jn| (EOS). Therefore, whether frame- 
dragging effects predicted by general relativity have ac- 
tually been observed in accreting neutron-stars depends 
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crucially on our knowledge of the moments of inertia of 
neutron stars. 

Theoretical calculations of neutron-star properties 
such as their masses, radii, and moments of inertia re- 
quire the knowledge of the equation of state of neutron- 
star matter up to densities 3> 10 15 gem -3 . However, 
the equation of state at densities higher than the nuclear 
saturation density (~ 2.7 x 10 14 gem -3 ) is still largely 
unknown and the subject of intense theoretical and ex- 
perimental research [[L3|. Therefore accurate predictions 
for neutron-star structural properties are hampered by 
the current uncertainty in the equation of state of high- 
density neutron-star matter flljj ]. 

In principle, the equation of state of matter at high 
densities could be constrained by astrophysical measure- 
ments of neutron-star masses and radii, although such 
measurements are relatively rare. The masses of the com- 
pact objects in double neutron-star systems measured via 
precise pulsar timing of the evolution of their binary or- 
bits 0] were found to cluster around ~ 1.35 Af© [ ^5[ . 
Mass estimates of neutron stars in other binary sys- 
tems depend strongly on the unknown inclination of the 
binary orbit and provide little additional information 
on neutron-star properties jl5[ ]. Astrophysical measure- 
ments of the radii of neutron stars based on the emit- 
ting area of their thermal emission during thermonuclear 
bursts [|16| , during the quiescence phase in transient sys- 
tems f |l T|| 7 or during the cooling phase of isolated neu- 
tron stars [[l8| are also difficult because of the systematic 
uncertainties in the predicted model spectra. Based on 
these measurements so far no significant constraints were 
imposed on any of the current equations of state. 

Bypassing the uncertainties of the equation of state at 
high densities, optimal bounds on the masses, radii, and 
moments of inertia of neutron stars have been derived 
under minimal assumptions about the validity of general 
relativity and the microscopic stability of neutron-star 
matter Stricter bounds can be obtained if the so- 

called causality limit is imposed, i.e., by requiring that 
the speed of sound is less than the speed of light every- 
where in the neutron star |2(|, although the validity of 
this requirement has been questioned. Such limits on the 
macroscopic properties of the neutron stars are largely 
independent of their unknown equation of state. 

In this paper, we obtain optimal bounds on the mo- 
ments of inertia of the neutron stars that show kHz QPOs 
under minimal or no assumptions for their structure and 
equation of state, taking into account the bounds on their 
masses and radii imposed by the observations of these 
QPOs. Our aim is to address the possibility that low- 
frequency QPOs in neutron star sources are related to 
general relativistic frame dragging which leads to Lense- 
Thirring precession of inclined orbits. 

In §2, we discuss our assumptions and method of 
solution of the relevant equations. In §3, we present 
numerical and semi-analytical bounds on neutron-star 
masses, radii, and moments of inertia and derive an EOS- 
independent maximum limit on the (Lense-Thirring) 



nodal precession frequency. In §4, we discuss the implica- 
tions of our results for the Lense-Thirring interpretation 
of the observed HBOs in accreting neutron stars. 

II. ASSUMPTIONS AND METHOD OF 
SOLUTION 

In the absence of detailed knowledge of the equation of 
state for matter at very high densities (~^> 10 14 g cm 3 ), 
we shall follow the procedure outlined by Sabbadini & 
Hartle to obtain bounds on the neutron-star masses, 
radii, and moments of inertia under a set of minimal as- 
sumptions for the equation of state of neutron-star matter 
above some fiducial energy density pq. These are: (a) the 
matter is cold, (b) the pressure in a given fluid element 
is determined uniquely by its energy density, (c) the en- 
ergy density and pressure are everywhere positive, and 
(d) dP/dp > (microscopic stability). In deriving these 
new bounds we make use of the recent observations of 
kHz QPOs in LMXBs to constrain the macroscopic prop- 
erties of neutron stars in these systems, i.e., their masses 
and radii ||]. 

The inferred spin frequencies of the neutron stars 
in the systems that show kHz QPOs and HBO are 
~ 250 — 350 Hz |||| , which are significantly smaller than 
the > 1500 Hz breakup frequency of a typical neutron 
star Jlij. We, therefore, assume that the neutron stars 
in all these systems are slowly rotating. Hereafter, we 
also set c = G = 1 , where c is the speed of light and G is 
the gravitational constant. 

The frequencies of the kHz QPOs observed in the X- 
ray brightness of accreting neutron stars and their de- 
pendence on mass-accretion rate have lead to the iden- 
tification of the highest-frequency QPOs with Keple- 
rian frequencies of stable circular orbits in the accretion 
disks ]3| H . In this interpretation, the observation of kHz 
QPOs from an accreting neutron star imposes two con- 
straints on its mass and radius (i) The radius of 
the orbit responsible for the highest-frequency kHz QPO 
must be larger than the radius of the innermost stable 
circular orbit around the neutron star. This leads to an 
upper bound on the neutron-star mass: 

M NS < (\/86W ma xr\ (1) 

where z^ max is the maximum observed frequency of the 
highest-frequency QPO. (ii) The radius of the orbit re- 
sponsible for the highest-frequency kHz QPO must be 
larger than the radius of the neutron-star itself. This 
leads to a mass-dependent upper limit on the neutron 
star radius: 

2M NS <i? NS <(Jg^Y /3 , (2) 

where the lower limit on the neutron star is simply the re- 
quirement that the central compact object is not a black 
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hole. In writing Eqs. (|l|) and (^) we have neglected the 
fact that the neutron star is slowly spinning; for the in- 
ferred spin frequencies of neutron stars in these systems, 
the correction is only < 10% towards increasing these 
upper bounds. 

In calculating the structure of a neutron-star with mass 
and radius consistent with the above bounds, we divide 
it into two regions, given a value of the fiducial den- 
sity po, above which we do not trust the equation of 
state: the core, with mass M c and radius R c in which 
/) > po, and the envelope exterior to the core. It has 
been shown |19|] that the combined stellar configuration 
satisfies the set of minimal assumptions for the equation 
of state for neutron-star matter if and only if 
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where Po is the pressure that corresponds to the energy 
density p at the edge of the core. 

Among all core configurations that satisfy the core 
mass limits (|^) and (Q), there is one that maximizes the 
stellar moment of inertia Jil| . The maximizing core con- 
figuration is that of constant density p c > pQ. For a given 
equation of state and fiducial energy density po , and for 
a given neutron-star mass and radius we calculate the 
structure of the neutron star for the core configurations 
that maximize its moment of inertia. We then scan the 
ranges of neutron-star masses and radii allowed by the 
kHz QPOs and obtain the global optimal bounds on the 
neutron-star moment of inertia as a function of po- 

To first order in the stellar angular velocity O, the 
metric in and around a slowly rotating star is |2jJ 



ds 2 = -e v ^dt 
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where u>(r) is the angular velocity of a locally non- 
rotating frame at radius r measured by an asymptotic 
inertial observer. Under these assumptions, the quanti- 
ties v(r) and m{r) satisfy the Oppenheimer-Volkoff equa- 
tions 
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as well as the equation 



dv ^ m + 47rr 3 P 
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where P is the pressure and p the energy density at a 
radius r inside the neutron star. In order to solve for the 
quantity co(r), and hence for the moment of inertia of the 
neutron star, we define, following Hartle 
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and 
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and solve the equation 
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dr 
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(11) 



for the radial dependence of /. 

Given a fiducial density po arL d the corresponding pres- 
sure Po, as well as a core mass M c and radius P c , we in- 
tergrate numerically Eqs. (|^) and (Q) outwards from the 
core-envelope interface until the intergration reaches the 
surface of the neutron star defined by P(Pns) = 0. The 
mass of the neutron star is then simply to(Pns) = Mns- 

For the core configurations that maximize the neutron- 
star moment of inertia (i.e., the one with constant en- 
ergy density in the core), we integrate analytically the 
same equations (Eqs. Q, 0) but inwards from the core- 
envelope interface. If the resulting neutron-star mass and 
radius satisfy the constraints imposed by the kHz QPOs, 
we then solve numerically Eq. ^) with the boundary con- 
dition i/(P NS ) = hi [1 - 2M ns /Pns]- 



Equation (11) is a second-order partial differential 
equation with the boundary conditions 



4f 

dr 







r=0 



and 



f(r = P NS ) = 1-2 



R 3 



(12) 



(13) 



where I is the moment of inertia of the neutron star given 
in closed form by j22\\ 
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dj(r) 
dr 



f(r)r 3 dr 
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'\- 1m(r)/r 



(14) 



We integrate Eq. (|11|) outwards from the center of the 
neutron star, using boundary condition ( p"2| ) and a trial 
value fo at the inner boundary. We call this trial solution 
/t r (r) and calculate the corresponding moment of inertia 
7 tr using Eq. (pr|). We note that Eq. ( |Tl| ) is scale-free 
and therefore assuming a different value of / at the inner 
boundary, f(r — 0) = £/o, where £ is a constant, results 
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in the solution f(r) = £/t r (0 and the corresponding mo- 
ment of inertia I = £/tr- Given the trial solution, we can 
then calculate the value of the parameter £ for which the 
boundary condition (113) is satisfied, i.e., 



Mr = Rns) + 



B 3 



It, 



(15) 



As a result, the solution of Eq. ( |Tl[ ) that satisfies bound- 
ary conditions (12) and (13) is just f(r) — £ftr(r) and 
the moment of inertia of the neutron star is 



I = 



Itr 



Mr = J?ns) + 27tr/i? 3 



(16) 



NS 



We performed all numerical integrations using a 
fourth-order Runge-Kutta scheme. We verified our nu- 
merical implementation of the procedure outlined above 
and our integration algorithms by taking the limit 
po — ► oo, i.e., assuming that we know the equation of 
state everywhere in the neutron star, and comparing the 
calculated moments of inertia for different neutron star 
parameters and equations of state with the values given 
by Cook, Shapiro, and Teukolsky Jl4]. We performed an 
additional test by taking the limit po — ► 0, i.e., assum- 
ing a constant-density neutron star, and comparing the 
calculated moments of inertia for different neutron star 
parameters with the values given by Sabbadini and Har- 
tle |L9[ and Abramowicz and Wagoner j2^]. In all tests 
the agreement was better than ~ 0.5%. 



III. BOUNDS ON NEUTRON-STAR 
PROPERTIES 

A. Optimal bounds on neutron-star masses and radii 

The constraints ([!]) and @ imposed on the neutron- 
star mass and radius by the identification of the highest- 
frequency kHz QPO with a Keplerian orbital frequency 
are derived without any assumptions regarding the equa- 
tion of state for neutron-star matter §. These con- 
straints can be further optimized assuming an equation of 
state for the envelope of the neutron star and constrain- 
ing the properties of the core (eqs. || and Q) so that 
they satisfy the minimal set of assumptions discussed 
in §2. For a given equation of state and fiducial den- 
sity po and for each core mass and radius in the allowed 
range defined by the limits @ and (0), we integrate the 
Oppenheimer-Volkoff equations and calculate the result- 
ing neutron-star mass and radius. 

Figure 1 shows the allowed ranges of stellar masses 
and radii for two representative equations of state (la- 
beled according to Ref. |Q) and for different values of 
the fiducial density pq. In order to understand the qual- 
itative behavior of the bounds shown in Figure 1, for a 
given equation of state, we call p c (M^s) the central den- 
sity and i? eos (AfNs) the radius of the star of mass Mns 
for this specific equation of state. 
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FIG. 1. Bounds on the mass and radius of a neutron star 
imposed by the identification of a 1220 Hz QPO with a Ke- 
plerian orbital frequency. Different panels correspond to dif- 
ferent equations of state, which were assumed to be valid up 
to some fiducial density. Different line types correspond to 
the constraints imposed for different values of the fiducial 
density. Dotted lines: equation-of-state independent limits 
(eqs. |lj and Q); short-dashed lines: po = 2.7 x 10 14 gem -3 ; 
long-dashed lines: po = 7 x 10 14 gem -3 ; solid lines: 
po = 2 x 10 15 gem -3 . The dotted lines for the maximum 
allowed radii in both panels overlap with the short-dashed 
lines and have been omitted for clarity. 



In the limiting case po — > 0, i.e., when the equa- 
tion of state is assumed to be unknown everywhere 
in the star, the lower bound on the neutron-star ra- 
dius corresponds to the general relativistic requirement 
-Rns > (9/4)Mns p9| , so that the central pressure in 
the neutron-star is not infinite. The upper bound on the 
neutron-star radius corresponds to the bound (^|) set by 
the kHz QPOs. 
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FIG. 2. Maximum neutron-star radius allowed by the iden- 
tification of a 1220 Hz QPO with a Keplerian orbital fre- 
quency, as a function of the fiducial density po up to which 
an equation of state was assumed to be known. Different 
curves correspond to different equations of state. The nuclear 
saturation density is p nB — 2.7 x 10 14 gcm~ 3 . 

For a non-zero fiducial density, the lower bound on 
the stellar radius corresponds to the requirement R c > 
(9/4) M c imposed on the core (and not the stellar) prop- 
erties. At high fiducial densities, the mass and radius 
of the envelope become non-negligible, i?Ns > -Rc, and 
hence the lower bound on the neutron-star radius for a 
given mass becomes tighter (moves to the right in Fig- 
ure 1). For a non-zero fiducial density, the upper bound 
on the stellar radius becomes tighter, as well. The ra- 
dius of a star of given mass decreases with increasing 
central density and hence, when po > Pc(^ns), the max- 
imum allowed radius of a neutron star with mass Mns is 
R eos (M-Ns), which may be smaller than the upper bound 
on the radius imposed by the kHz QPOs. 

In Figure 2 we plot the maximum allowed neutron- 
star radius as a function of the fiducial density, for dif- 
ferent equations of state (labels according to Ref. [fl4j|). 
As indicated from Figure 1 already, using any equa- 
tion of state up to about the nuclear saturation den- 
sity (p ns ~ 2.7 x 10 14 gem -3 ) does not not alter the 
bounds on the neutron-star radius (eq. |^]) imposed by 
the kHz QPOs. Extending the use of any equation of 
state to higher densities (up to two or four times p ns ) 
leads to tighter constraints on the stellar radius that de- 
pend strongly on the assumed equation of state. The 
bound imposed by Eq. (|l|) on the neutron star mass is 
not affected by the knowledge of the equation of state. 

Note here that all physical equations of state satisfy 
the set of minimal assumptions used in constructing Fig- 
ures 1 and 2. As a result, the bounds on the mass and 
radius of a neutron star plotted in these figures will not 
lead to any constraints on the microscopic properties of 
the equation of state at high densities in addition to those 
already imposed by relations M) and (0) . However, these 



considerations do impose additional constraints on the 
macroscopic properties of individual neutron-stars, which 
in principle can be compared with other independent es- 
timates of their masses and radii. Such additional esti- 
mates may be obtained from the X-ray spectra |l6| ] or 
the amplitudes of the coherent oscillations during X-ray 
bursts Q. 

B. Bounds on neutron-star moments of inertia 

For a given equation of state and fiducial density po 
and for a neutron star with mass and radius consis- 
tent with the bounds calculated above, we calculate the 
upper bound on the stellar moment of inertia, as out- 
lined in §2. Figure 3 shows the maximum value of the 
quantity I i5 /(M^s/M Q ), where / = Z45IO 45 gem 3 is 
the neutron-star moment of inertia, as a function of the 
fiducial density for different equations of state and for 
a maximum kHz QPO frequency of 1220 Hz. The maxi- 
mum ratio I^/{M^s/Mq) decreases with increasing fidu- 
cial density po, in part because //Mns -Rns anc ^ the 
maximum allowed neutron-star radius decreases with in- 
creasing po (see Figure 2). For fiducial densities higher 
than the nuclear saturation density, the maximum ra- 
tio I^5/(Mns/Mq) depends on the assumed equation of 
state. 

For a neutron-star showing kHz QPOs with maximum 
frequency i/ max , we can obtain the maximum value of 
the ratio I45 / (M^s / Mq) , independent of the equation of 
state, by taking the limit po — > 0. This corresponds to a 
star with constant density and is the configuration that 




FIG. 3. Maximum neutron-star moment of inertia in units 
of its mass allowed by the identification of a 1220 Hz QPO 
with a Keplerian orbital frequency, as a function of the fidu- 
cial density po up to which the equation of state was as- 
sumed to be known. Different curves correspond to differ- 
ent equations of state. The nuclear saturation density is 
p ns = 2.7 x 10 14 gem -3 . 
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FIG. 4. Maximum neutron-star moment of inertia in units 
of its mass as a function of the maximum observed kHz QPO 
frequency. The upper bounds calculated using the numerical 
method described in §2 (solid line) are compared to the ones 
calculated analytically in the limit of zero (dashed line) and 
maximum (dotted line) compactness of the neutron star. 

maximizes its moment of inertia. Figure 4 (solid line) 
shows the resulting dependence; for a maximum kHz 
QPO frequency of 1220 Hz, I 45 /(M NS /M G ) < 2.3. 

We can obtain the dependence of the maximum ratio 
I^/(M^/Mq) on fmax a t the limit po — > in the fol- 
lowing analytical way. For any star, independent of the 
equation of state, Sabbadini and Hartle fill showed that 



/ < £M NS R 2 NS 



(17) 



The parameter £ is a function of the compactness of the 
star (Mns/-Rns) an d has the limiting values £ m i n = 2/5 
in the Newtonian limit (Mns/-Rns ~~ * 0) and £ max = 
0.799 ~ 4/5 in the limit of maximum compactness al- 
lowed by general relativity (i.e., Mns/-Rns = 4/9) jl9). 
Combining equation (|l7]) with the radius bounds (Q) and 
(0), we obtain 



A I] 



< 



< 



i 



NS 



24^i ax " 307r2^ ax 



(18) 



where £ corresponds to the compactness of the neutron 
star that maximizes the ratio I /Mrs and in the last in- 
equality we used the fact that £ < 4 /5 . In Figure 4 we 
compare the analytical upper bound (|l8| ) for the two lim- 
iting values of £ (£ = 2/5: dashed line; £ = 4/5: dotted 
line) to the one calculated numerically in the limit po — > 
(solid line). 

When po — > 0, the maximum value of the ratio 7/A/ns 
corresponds to a neutron star of constant density with 
typically the maximum radius allowed by constraint (|^) 
imposed by the kHz QPOs. For such a neutron star, 
M NS /i? NS = 1/6 (see Eqs. jjj and (|]) and hence £ is 
closer to the Newtonian limit of 2/5 than to its maxi- 
mum value of 4/5, as suggested by Figure 4. Indeed, by 



comparing the analytical scaling to the numerical result 
we find that a better approximation to the upper bound 
of the ratio I45 / (Mhs/Mq) corresponds to £ = 0.452 and 
is (see also Figure 4) 



J 45 
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(19) 



This is the maximum value of the ratio '(Mns/M©) 
allowed by general relativity for a neutron star that shows 
a maximum kHz QPO frequency of v max in its power 
spectrum. 



C. Maximum nodal precession frequency 

The orbital plane of an infinitesimally inclined circu- 
lar orbit around a slowly spinning neutron star precesses 
because of general relativistic frame dragging, as well as 
because of classical effects related to the quadrupole mo- 
ment of the stellar gravitation field. The nodal precession 
frequency j^np of such an orbit is ]10[ 



NS 



(20) 



where z/lt and vq are the contributions of the general 
relativistic Lense-Thirring precession and of the classical 
precession, respectively, v s is the spin frequency of the 
neutron star, and vk is the Keplerian orbital frequency 
of the orbit. 

The frequency of the general relativistic precession is 
directly proportional to the ratio I/M, for which we 
have obtained optimal bounds in the previous section. 
Given that the effect of classical precession is to re- 
duce the nodal precession frequency and using the upper 
bound (|l8|) on the ratio I/M^s we obtain 



45.2 



300 Hz 



Hz , 



(21) 



where we used £ ~ 0.452, for the reasons discussed in 
§3.2. 

Equation ( ^l| ) shows that there exists an upper limit on 
the nodal precession frequency of a circular orbit around 
a slowly spinning neutron star, which we have obtained 
analytically. It is remarkable that this upper limit is in- 
dependent of all the other properties of the neutron star, 
of the unknown equation of state, or of the properties of 
the circular orbit. 



IV. DISCUSSION 

We have calculated optimal bounds on the masses, 
radii, and moments of inertia for slowly rotating neutron 
star in which observed frequencies of kHz QPOs have 
been identified with Keplerian orbital frequencies. As- 
suming the validity of an equation of state up to some 



G 



fiducial density po , these bounds become tighter as po in- 
creases. In the limiting case of po — > 0, i.e., when we make 
no assumption regarding the equation of state, we have 
derived an analytical upper bound on the neutron-star 
moment of inertia for a given maximum observed kHz 
QPO frequency. We also obtain analytically the maxi- 
mum nodal precession frequency of an inclined circular 
orbit around a neutron star, which depends only on the 
spin frequency of the star and is independent of the other 
stellar properties, the equation of state, or the properties 
of the circular orbit. 

In this section we use these constraints to address 
the possible observational evidence for general relativistic 
frame-dragging effects in the rapid variability of accret- 
ing neutron stars (h]]. Many such sources often show 
three distinct types of QPOs that are not harmonically 
related The two QPOs at kilohertz frequencies are 

believed to occur at the Keplerian frequency of a stable 
circular orbit in the accretion disk and at its beat with 
the neutron-star spin frequency j^H . According to Stella 
& Vietri Jh}, the third, low-frequency QPO (the HBO) 
occurs at the nodal precession frequency of the orbit re- 
sponsible for the kHz QPO. 

We can first test quantitatively the suggestion that 
the HBO occurs at the nodal precession frequency of 
an inclined circular orbit by comparing the maximum 
observed HBO frequency in different sources with the 
maximum possible nodal precession frequency around a 
slowly spinning neutron star (Eq. pl|). We will use the 
data of five bright neutron stars, in which all three QPOs 
can be identified unambiguously, as discussed in detail in 
Ref. pHl . Figure 5 shows the maximum HBO frequencies 
in these sources plotted against their spin frequencies in- 
ferred from the peak separation of the kHz QPOs. The 
data points are compared to the maximum possible nodal 
precession frequency around a neutron star calculated in 
§3.2 (Eq. (2^]). Four out of the five sources are incon- 
sistent with our optimal bound at least at the one-sigma 
level. In particular, the observed HBO frequencies in the 
source GX 17+2 can be excluded at high statistical sig- 
nificance from being nodal precession frequencies, if the 
peak separation of the kHz QPOs is equal to the spin 
frequency of the neutron star. 

It is possible that the observed HBO occurs at the 
second harmonic of the nodal precession frequency. A 
precessing circular orbit has a two-fold symmetry that 
could, in principle, produce even-order harmonics that 
are stronger than the odd-order harmonics. The observa- 
tion of a subharmonic of the HBO in neutron-star sources 
as well as of subharmonics of similar QPOs in black-hole 
sources |25| gives additional weight to this conjecture. In 
this case, we address the nodal-precession interpretation 
of the HBO using the observed correlation between the 
HBO and kHz QPO frequencies in the same five bright 
neutron-star sources discussed in Ref. JlTj ] . 

When the effects of classical precession are negligible, 
observation of a Keplerian, spin, and nodal precession 
frequencies leads to a direct measurement of the ratio 
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FIG. 5. Maximum HBO frequency versus inferred spin fre- 
quency for five bright neutron-star sources. The dashed line 
is the maximum nodal precession frequency for an inclined 
circular orbit around a star of a given spin frequency. 

I/M of the neutron star (see Eq. ^0|). Identification 
of the HBO frequency in four of these five sources with 
the second harmonic of the nodal precession frequency 
requires I45 / (M^s / Mq) > 2.3 at the 99% confidence 
level [[n]. However, for the maximum observed kHz QPO 
frequency of 1220 Hz, I 45 /(M ns /Mq) < 2.3, indepen- 
dent of the equation of state (extreme case of po — 0). 
Moreover, assuming even the unrealistically stiff equa- 
tion of state L up to twice nuclear saturation density 
results in an upper limit of J45 / (Mns / Mq) < 1.7 (see 
Figure 3). Assuming any other equation of state results 
in even lower values of the ratio I45 / (Mrs / M@) which 
are inconsistent at high statistical significance with the 
value required by the nodal-precession interpretation of 
the HBO. 

We therefore conclude that the nodal-precession in- 
terpretation of the HBO observed in several neutron- 
star sources is inconsistent with the identification of the 
higher-frequency kHz QPO with a Keplerian frequency 
of a circular orbit and with the identification of the fre- 
quency separation of kHz QPOs with the spin frequency 
of the neutron star ETJl . 
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